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Abstract 

Exactly solvable rererence potentials of several smoothly joined Morse-type components were 
constructed for the lowest two excimer states of Ai2 molecule. The parameters of the potentials 
have been ascertained by fitting to the experimental data, and they are reliable in a wide range 
of nuclear separations (r > 1.9 A). A large number of quantum mechanical Franck-Condon factors 
for the 0+ — ► 0^" and l u — > 0+ transitions have been calculated and compared with the observed 
spectroscopic features. The fitting procedure also involved dipole transition moments, which have 
been adjusted to the known radiative lifetimes of the vibrational levels. The resulting potential 
energy curves accurately reproduce the first and the second emission continua of Ar^j as well as the 
oscillatory spectrum related to their inner turning point region. The numbering and the positions 
of the vibrational levels reported by Herman et al. [P. R. Herman, P. E. LaRocque, and B. P. 
Stoicheff, J. Chem. Phys. 89, 4535 (1988)] have been confirmed. The calculated spectroscopic 
parameters, w e = 296.26 cm^ 1 , D e = 6128.3 cm -1 for the 0+ state, and w e = 287.30 cm , 
D e = 5929.6 cnr 1 for the l u state, are consistent with previous estimations. Separation between 
the minima of the potential wells at R e = 2.3893 A was found to be 50.60 meV, compared with 
75.24 meV splitting in the separated-atom limit. 
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I. INTRODUCTION 



Rare gases (RG) are important laser media in the vacuum ultraviolet (VUV) spectral 
region [J Q] , and this is one of the reasons for the undying interest in these chemically 
inactive species. Normally, they are almost atomic gases, because electronically unexcited 
RG dimers are only weakly bound (except radioactive Rn2 with dissociation energy of about 
400 K {3I). In contrast to this, the electronically excited RG dimers (excimers) in their 
lowest states l u and 0+ are strongly bound. These excimers can be created with the help of 
various excitation sources (charged particle beams, synchrotron radiation, etc.), and their 
decay produces photons in the VUV region, forming two intense continua observed long ago 



for all RGs (except Rn) Q, J E, (see also 0, [u 



for the recent Ne anc. 
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A vast literature exists on the RG dimers in their ground electronic state 
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1 171 Il8| ] , because they are good prototypes for studying van der Waals interactions. 
At the same time, there is still a lack of knowledge about the RG excimer states, which 
are of crucial importance for the development of efficient light sources in the VUV region 
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231 ] . Of special interest in this context are the heavier excimers for 
which the laser emission has been reported The observed laser effect is related to 

the so-called second continua, broad structureless emission bands with maxima at 126, 146 
and 172 nm for Ar^, Kr^, and X.e^, respectively. They correspond to transitions from the 
vibrationally relaxed l u and 0+ states to the repulsive part of the ground electronic state 0^", 
in contrast to the so-called first continua originating from the highest vibrational levels of 
the 0+ excimer state, and therefore being spectrally close to the atomic 3 Pi — > 1 S resonance 
lines. The first continuum is related to the excimers' classical outer turning point region, 
while the weak oscillatory emission from the region of their inner turning points has also 
been observed under selective synchrotron radiation excitation 24]. 



All these important spectroscopic features can be uniformly treated [25|, [26j, but to this 
end, one needs reliable potential energy curves (PECs) and coordinate dependencies of the 
related transition moments, as well as a reasonable model to describe vibrational relaxation 
of the excimers. In this paper, we concentrate on the first two issues, while the relaxation 
dynamics and the time- dependent emission spectra will be studied in a separate paper 27^ . 
Valuable information on the energetic position, absolute numbering and radiative lifetimes 
of the higher vibrational levels for all heavier RG excimers has been obtained in a series 
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of supersonic jet expansion experiments 28, 
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30 



3JJ , using tunable VUV radiation to 



produce high-resolution fluorescence excitation spectra. Unfortunately, these high-quality 
reference data concern only a narrow range of vibrational levels near the dissociation limit, 



and one cannot construct reliable PECs for RG excimers sole 
and semi-empirical potentials are available 
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y on this basis. Various ab initio 
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371 . |38| , which satisfactorily 



reproduce some of the characteristic spectroscopic features, but not the whole variety of 
available experimental data. Apart from imperfections of the PECs of the excimer states, the 
discrepancies between theoretical and experimental spectra may be caused by inaccuracies 
of the ground state potentials in the short distance range. 

Thus, the spectroscopic properties of RG dimers certainly need further investigation. 
This is the motivation for the present study, the main goal of which is to construct accurate 
reference potentials of relatively simple analytic form, by directly fitting their parameters to 
the observed data, and taking account of the relevant theoretical considerations. In principle, 
the approach that will be described below can be applied to any RG dimers, since they have 
much in common. For several reasons, however, only Ar2 molecules are under examination in 
this paper. First, compared with Xe and Kr, the effects of spin-orbit coupling are much less 
pronounced in Ar, which makes the theoretical analysis more simple and, hopefully, more 
adequate. Second, a sufficient amount of good-quality experimental data are available for 
Ar gas, which, again, is an advantage compared with Xe and Kr. Third, there is a practical 
point of special interest in Ar^ excimers, as they produce the highest energy photons (9.8 
eV) for the potential laser applications, which is the field of great interest and permanent 
development 0, I40}. It therefore seems reasonable to test the method on Ar2 molecules, 
leaving more complex cases for the future analysis. 

The reference potential approach to different spectral problems has been described in 
detail elsewhere (see, e.g., 25|, I4JJ and references therein). In this paper, we are going 
to construct smooth multi-component exactly solvable reference potentials for the lowest 
excimer states l u and 0+, and for the ground electronic state 0+ of the Ar2 molecule. These 
approximate PECs should be realistic in a wide range of nuclear separations related to the 
observed spectroscopic features. To achieve this goal, fully quantum mechanical calculations 
of bound-bound and bound-free Franck-Condon factors for l u — > 0+ and 0+ — > 0+ transitions 
have been performed. The procedure has to be repeated many times, varying the parameters 
of the PECs and the transition moments. In this way, step by step, the calculated energy 
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levels and their radiative lifetimes, as well as the intensity patterns of Franck-Condon factors, 
can be fitted to the available experimental data. 

The paper is organized as follows. In Section II, a short overview of the theoretical 
background of the method is given, specific properties of the reference potentials for the 0+, 
0+ and l u electronic states of Ar 2 are described (in sub-sections A, B, and C, respectively), 
and the details of the fitting procedure are explained. Section III is devoted to demonstrating 
the calculated Franck-Condon factors, transition probabilities and radiative lifetimes for the 
vibrational levels of the 0+ and l u states. Finally, a brief discussion of the obtained results 
and future prospects of the research are given in Section IV, which concludes the paper. 



II. EXACTLY SOLVABLE MULTI-COMPONENT REFERENCE POTENTIALS 



At first sight, the reference potential approach we are going to adopt may seem naive and 
unjustified. However, the solution to a problem always depends on the context. Modern 
spectroscopic analysis often involves many thousand items of high resolution input data, 
but unfortunately, such detailed information is not yet available neither for the ground 
state of Ar 2 molecule nor for the excimer states l u and 0+. Indeed, the scattering data for 
these species are almost lacking. Moreover, even the total number of the bound states is 



not conclusively established, and on 
ascertained for the excimer states 



y a narrow range of discrete levels has been accurately 



301 ]. while their absolute numbering still needs to be 



confirmed. Under these circumstances, there is no chance to deduce the "real" PECs directly 
from the experimental data. One can only take some steps towards this goal by combining 
the available experimental data with the relevant theoretical results. In this context, the 
simple but mathematically rigorous approach in question has proven to be quite useful. 

To begin with the analysis, let us recall another important peculiarity of the system 
under study. Namely, for all RGs, there is a big difference (about 1.4 A in the case of Ar 2 ) 
between the minima of the tiny potential well of the ground state and those for the excimer 
states. It means that bound-free transitions to the repulsive wall of the ground state PEC 
are responsible for the major part of the observed fluorescence spectra. Of course, as we 
demonstrate, bound-bound transitions are also very important to reveal the details of the 
first emission continuum. Unfortunately, this fine structure, which should contain tens of 
spectral lines in the range of less than 0.1 eV width (compared with more than 4.5 eV width 
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of the overall spectrum), has not yet been detected experimentally. Moreover, the spectral 
resolution in fluorescence experiments on Ar and other RGs is typically ~10 meV. Therefore, 
at the present level of our knowledge, one should not put too much effort into deducing the 
PECs that would reproduce the observed (relatively few) level positions with the utmost 
accuracy. Formally, one can always achieve this specific goal, but this does not guarantee 
the reliability of the potentials in the major part of the domain. Indeed, as is well known, 
even the full knowledge of the discrete energy spectrum combined with the full scattering 



information is insufficient to uniquely determine the potential [421 ]. 

In view of the above, we set the following criteria to the quality of the reference poten- 
tials for the Ar?; excimer states: 1) these PECs should reproduce the observed bound-bound 



transition energies 30| with ~1 meV accuracy; 2) they should reproduce all details of the 
fluorescence spectra in the range from about 7 to 11.63 eV with ~10 meV accuracy. To 
achieve these goals, we use a fitting method which is based on repeated and accurate calcu- 
lation of the Franck-Condon factors. Naturally, this strictly quantum mechanical procedure 
presumes highly efficient solution of the Schrodinger equation, because this elementary act 
has to be performed many million times within a reasonable time scale. For this reason, we 
try to adopt an analytic approach of calculating the energy eigenf unctions, which indirectly 
means that we have to construct exactly solvable approximants (reference potentials) to the 
"real" PECs. If one sets a "physical constraint", requiring continuity of the potential and 
its first derivative in the whole physical domain, the number of suitable options becomes 
quite limited. Certainly one of the best choices is to compose a reference potential of several 
Morse-type [43( pieces 

V(r) = V k + D k [exp(-a fc (r - R k )) - l] 2 , re (0, oo), (1) 

where V k , D k , a k and R k are real (not definitely positive!) constants, and the subscript 
k = 0, 1, 2. ..labels different components smoothly joined at the boundary points X k+ i. 

As is well known, the classical Morse potential belongs to the family of shape invariant 
potentials 44|, and its energy eigenvalue problem can be solved with the help of solely 
algebraic techniques. The shape invariance is lost, if there are several analytically different 
components, as assumed in Eq. (1). Consequently, the discrete energy levels of such multi- 
component potentials cannot be given in an explicit analytic form. With some concession, we 
can still preserve the term "exactly solvable" for this kind of piece-wise potentials, because 
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the two linearly independent solutions of the related Schrodinger equation can be always 
found analytically to any desired accuracy. Indeed, introducing a dimensionless variable 
y k = 2a k exp(— a k {r — Rk)), the Schrodinger equation for a Morse-type PEC can be converted 



into a confluent hypergeometric form 45| 



d G(a k , fi k ;y k ) . /„ . ->dG(a k , fi k ;y k ) , 
Vk -t~2 H ( 2 A*fc + 1 - Vk) j H K - /xjb - l/2)Gr(a fc , yk) = 0. (2) 



Here a k = \J D k /C /a k} C = j- is a characteristic constant (C = 0.10464 meV-A 2 for Ar2), 
fi k = ^/ <yVk+Dk E ^ c ; anc i the solutions of the Schrodinger equation read (up to normalization) 



$ = exp(—y k /2)y k k G(a k , fi k ; y k ). Note that the parameters a k , \i k and the coordinate y k may 
be imaginary, if D k < or the total energy E > V k + D k . 

There are several possibilities to construct the fundamental system of solutions for Eq. 



(2) (see 45[], Chapter 6, for a thorough overview). For example, one can make use of the 



special solutions 

G?! = $(-a* + + 1/2, 2fi k + 1; y fe ), G 2 = y~ 2Atfc $(-a fc - /i fe + 1/2, -2// fc + 1; y fc ), (3) 
where the symbols 

_ , . 1 a - x a(a + 1) ■ x 2 . A . 

denote confluent hypergeometric functions. Correspondingly, the two linearly independent 
solutions of the original Schrodinger equation read 

*ifci = y k h S(a k ,fi k ; y k ), * k2 = y k ^ k S(a k , -fi k ;y k ), (5) 
where we have introduced another very useful function 

S(a, c; x) = exp(-x/2)$(-a + c + 1/2, 2c+ 1; x), (6) 
which can be evaluated as follows 461 ]: 



S(a,c;x) = X>„, B = 1, B l = -^y, ( 7 ) 

n=0 

S " = /o X , — r(-a5„,-i + 7^-2), n = 2, 3, ... 
n(2c + nj 4 

From Eqs. (5) to (7) one can infer that ^ k i and ^ k 2 are always complex conjugates, if 
E > V k + D fc . 
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A. Four-component reference potential for the ground state of Ar2 



After the brief theoretical introduction, let us try to put the ideas into practice. Our 
fitting procedure does not actually involve the parameters of the ground state, but only 
those of the excimer states l u and 0+. In other words, the PEC for the ground state is 
assumed to be fixed. Nevertheless, in order to solve the whole quantum mechanical problem 
exactly, we have to construct an exactly solvable (in the above-mentioned sense) reference 
potential for the ground state as well. Of course, such a constraint is technical rather than 
conceptual. The Schrodinger equation can be always solved numerically, but in the present 
context the main problem is how to fix a sufficiently realistic PEC. The shape of the ground 
state potential for Ar 2 is accurately known in the intermediate and long-distance range, but 
the knowledge about its repulsive short-distance part is rather ambiguous. As mentioned, 
this region is very important to interpret the observed spectroscopic features, thus our aim 
is to describe the repulsive wall as adequately as possible. 

The resulting reference PEC for Ar 2 is shown in Fig. 1, and its parameters are given in 
Table 1. All four components have the simple analytic form of Eq. (1), and their parameters 
have been ascertained from the least squares fit to the ab initio CCSD(T) daug-cc-pV5Z- 
33211 potential by Fernandez and Koch 47|. This particular PEC has been preferred, 
because this is the only work to date where both the ground state and the excimer state (0+) 
potentials have been calculated on the same theoretical basis. These PECs nicely reproduce 
the observed spectra related to the excimers' inner turning point region [2^, which is a 
strong argument for their reliability at short distances. As can be seen in Figs. 1 and 2, the 
ground state PEC by Fernandez and Koch is in good agreement with the other available ab 
initio potentials 18|, |48( . The reference potential for Ar 2 is almost indistinguishable from 
the original PEC on the scale used in these figures. Therefore, it is expected to be reliable 
enough for our purposes. 

Now, let us have a look at the structure of this reference PEC. The two most internal 
components (r < Aj an d r e [Ai,X 2 ], respectively) represent the so-called pseudo- Morse 
(PM) potentials 25|, |4l| smoothly joined at the boundary point X\. Their important pecu- 
liarity is that the parameters Dk and a& {k = 0, 1) are not independent, but = \Ca\. 
It means that the tiny potential well is just of the limit depth to entirely lose the discrete 
energy spectrum. Since = 1/2, and the pseudo- Morse approximation is used only in the 



7 



region where E > V k + D k , the two special solutions ^ki and ^ k2 of the Schrodinger equation 
are complex conjugates. According to Eqs. (5) to (7), 

v£r fcl = yf"S(l/2,ip k ;y k ) = A k {y k )e iBk ^ e'^^, (8) 

where 

A k (y k )e iB ^ = 1 - Q VklA , (9) 

yy ' if3 k + 1/2 v ; 

■ (y^f L _ _vg±\ (^/4) 4 / _ y k /A \ 

(ifa + 1/2) 1! V i/3 k + 3/2 J (i/3 h + 1/2) (ift + 3/2) 2! V + 5/2 J 

and, consequently, the general solution reads 

*fc(r) = N k A k {y k ) cos [Bfc(j/ fc ) + <^ fc - a fc /5 fc r] , k = 0, 1. (10) 

Here = |/ttfc| , while the normalization factor A& and the phase constant (p k should be 
determined from the continuity requirements of the wave function and its derivative. 

The phase constant ipo for the most internal PM component can be always (in most cases 



quite easily) ascertained analytically 41], taking account of the physical boundary condition 



^ — ► as r — > 0. The next phase constant (pi can then be determined from the boundary 
condition = ■ In all cases analyzed in this paper the wave function's logarithmic 

derivative can be expressed in the form 

Sfc(r) tany? fc + T k (r) 



R k (r) + 



(11) 



U k (r) tcLiiip k + W k (r)_ 

where R k , S k ,T k ,U k and W k are some characteristic functions that do not depend on phase 
constants. Thus, if y?o is known, (pi can be easily determined. 

The region r G [X 2 ,X3] (including the minimum of the PEC) is approximated by an 
ordinary Morse (OM) potential (k = 2). Again, since the energy range E < V 2 + D 2 is 
out of interest, the special solutions \l/ 2 i and \l/ 22 are complex conjugates, and the general 
solution of the Schrodinger equation becomes 



where 



Mr) = N2[C 2 {r)tfm<p 2 + D 2 {r)], (12) 

C 2 (r) = Re [S(a 2 , i(3 2 ; y 2 ] sm(k 2 r) - Im [S(a 2 , i/3 2 ; y 2 ] cos(/c 2 r), 
D 2 {r) = Re[S(a 2 ,i/3 2 ;y 2 ] cos(k 2 r) + Im [S(a 2 , i(3 2 ; y 2 ) sin(/c 2 r), 
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k 2 = a 2 /? 2 , and fa = l/^l • To ascertain the phase constant <f2, one uses Eq. (11) and applies 
the boundary condition §^g} = §f{§^. 

The most external region r > X3 (A; = 3) is approximated by a reversed Morse (RM) 
potential with a negative "dissociation energy" D 3 . This might seem unphysical, because 
there is actually no hump on the original PEC. However, since the hump of the reference 
PEC is very small and located in the long-distance range (see Table 1) where V(r) pa 0, 
this artificial effect is nearly negligible for spectroscopic applications, while the analytic 
treatment remains as simple as in previously studied cases. Of course, one can introduce 
more such components and gradually shift the maximum to an arbitrarily long distance (in 
principle, to infinity), thus practically eliminating the artificial barrier. In this paper, this 
physically motivated but tedious procedure has not been undertaken, because the desired 
quality of the PECs (see the criteria set in Section II) can already be achieved with the help 
of only 3-4 components. 

A thorough analysis of the bound-states region E < V 3 + D 3 = has been given elsewhere 



49j . For the scattering states (E > 0), the wave function reads 



2C 3 (r) 



* 3 (r) = ={cos [D 3 (r) - kr] - tan ^3 • sin [D 3 (r) - Ax]}, (13) 

A/tam (p 3 + 1 



where C 3 (r) exp [iD 3 (r)] = S(ia 3j ifa;ix 3 ), f3 3 = \n 3 \,x 3 = \y 3 \ , and k = a 3 (3 3 = y/E/C. 
Since C 3 (r) — > 1 and D 3 (r) — > as r — > 00, ^ 3 (r) asymptotically approaches the free- 
wave form, ^(r) pa 2cos(y93 — kr). Consequently, the main spectral characteristic of the 
scattering states, the phase shift, reads S(k) = (n + l/2)7r — ip 3 , where n is an integer. 
To ensure the correct energy normalization, Eq. (13) has to be multiplied by the factor 

/ \-i/2 rn 

F = ( Am\J EC J 50|. As previously, to ascertain the phase constant ip 3 , one has to 
suitably adjust Eq. (11) and apply the boundary condition §^73^ = ^lexl ■ Thereafter, 
one can fix the normalization factors A^, N\ and Ao, using the continuity conditions for the 
components of the wave function 

^ 3 (X 3 ) = ^ 2 (X 3 ),^ 2 (X 2 ) = V 1 (X 2 ),* 1 (X 1 ) = * (^). 

Thus we have explained all details of calculating the stationary wave functions for an 
exactly solvable multi-component reference potential. The described scheme remains the 
same, independent of how many analytically different components one introduces. In view 
of the incomplete knowledge about the real PEC for Ar 2 , including just four smoothly joined 
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Morse-type pieces seems quite optimal. In fact, the full discrete energy spectrum E n < 
(n = 0-r6) can be accurately ascertained with the help of only two components (OM + RM), 
because the region r < X 2 = 3.35 A becomes practically negligible for these calculations. 
Determining the energy levels therefore reduces to a very simple zero-finding problem for a 
function, which is uniquely determined by the physical boundary conditions ^2{E,r) — > 
as r — > and ^s(E, r) — > as r — > oo, complemented with the condition ^^'^ = §4^] for 
an arbitrary r > X 2 (e.g., r = X 3 ). These three conditions can be fulfilled simultaneously 
only for the true energy eigenvalues, i.e., if E = E n . 

In Fig. 1 one can see little discrepancies between the original PEC and its approximant, 
which could be easily reduced by adding more components. An important point to discuss 
in this context is the physically correct long-range behavior of the potential. Indeed, apart 
from the artificial potential barrier, RM approximation seems to be absolutely incompatible 
with the attractive inverse power series coordinate dependence, which is expected at long 
distances. However, in the framework of the proposed approach, the "unphysical" nature 
of the RM approximation and the discrepancies mentioned are not substantial. First, as 
can be seen in Table 1, the calculated vibrational levels for the reference PEC fit with the 
observed ones 



30, 



51] even better than those of the original potential. Thus, adding more 



components to the potential is not motivated. Second, in a wide range of actual interest the 
RM approximation does not contradict to the inverse power series expansion (see the inset 
of Fig. 1). Therefore, since the described approach is reliable enough, and provides simple 
analytic solution to the problem, it has been used in the long distance region as well. 

As mentioned, the parameters of the ground state reference potential were kept un- 
changed. A large number of scattering (E > 0) wave functions for this PEC have been 
calculated, along with the wave functions of the bound states. These eigenfunctions are 
needed as a basis for further calculations, and this basis should be sufficiently complete to 
accurately reveal all details of the Franck-Condon factors for the vibrational levels of the 
excimer states. 



B. Reference potential for the 0+ state of Ar^ 

Construction of the reference potentials for the excimer states can be performed in the 
manner described in the previous subsection. The main difference is that the parameters 
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of the components are now treated as variables to be fitted to the experimental data. In 
addition, the ab initio CCSR(T) daug-cc-pV5Z potential for the 0+ state [4^] in the range 
r < 2.7 A has also been used as an input for the fitting procedure. Therefore, in this region 
(including the minimum at 2.3893 A) the reference potential practically does not differ from 
the PEC calculated by Fernandez and Koch. 

As previously, the reference potential was built up of several smoothly joined Morse-type 
pieces. Since spectroscopic applications are related to the excimers' bound states, the short- 
distance region (r < 2 A) is of less interest, and therefore, the reference PEC constructed 
for the 0„ state does not include any PM components. Thus, the resulting curve shown in 
Figs. 2 and 3 has only three constituents with parameters given in Table 2. The region 
r < Xi was approximated by an OM potential, and two RM components (r G [Xi, X 2 ] and 
r > X2, respectively), have been introduced for the long-distance range where spin-orbit 
coupling becomes important. Note that we cannot use the ab initio PEC by Fernandez and 
Koch for this range, because spin-orbit coupling was ignored in their calculations. 

As mentioned, the PM component remained nearly unchanged, while the RM compo- 
nents have been largely varied, preserving the continuity of the potential and its derivative, 
and trying to achieve a good fit with the experimental data. First, for any intermediate ref- 
erence PEC, the energy eigenvalue problem has been solved. As explained in the end of the 
previous sub-section, the discrete energy levels can be always found as the zeros of a charac- 
teristic function, which is uniquely determined by the physical boundary conditions and the 
continuity requirements (see, e.g., 49| and references therein for more details). Thereafter, 
the Franck-Condon factors for all vibrational levels were calculated and integrated over the 
full energy range. The results have been compared with the corresponding experimental 
data. If needed, the parameters of the reference PEC have been slightly changed and the 
whole procedure has been repeated until the desired quality of the fit was achieved. 

Naturally, the correct PEC should reproduce the second continuum and the real disso- 
ciation limit (E a = 11.623592 eV for the atomic 3 Pi level 52]). In addition, the fitting 
procedure involved the observed level positions (30I and the intensity patterns from the 
excimers' inner turning point region [24]. Unfortunately, only rough estimations for the 
important spectroscopic constant u) e of the 0+ state are available (see, e. g.. J30lp . but it has 
to be close to the values reported for the l u state {ui e = 293 ± 4 cm -1 [53j, oj e = 299 ± 3 



cm 1 



54j) and for the 2 T^, 2u state of Ar^ {oj e = 307 ± 0.4 cm 1 |55j). Finally, one has to 
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take account of the known radiative lifetime for the 0+ constituent of the second continuum, 



t = 4.20 ±0.13 ns 



561 ] . To achieve agreement with this value, the transition moment given 



by Fernandez and Koch has been corrected in the range r < 3 A (see the inset in Fig. 2). 
Table 2 and Figs. 3 and 4 demonstrate the obvious success of the described procedure. 



C. Reference potential for the l u state and l u — ► 0+ transition moment 

The reference potential constructed for the l u state looks similar to the 0+ PEC just 
described (see Fig. 2), although the fitting procedure was slightly different. This time, no 
ah initio PECs can be used for comparison, but on the other hand, a lot more experimental 
data are available. In addition to the level positions and their absolute numbering in the 
range v ' — 23 -r- 31 0, the radiative lifetimes for v' — (r = 3.2 ± 0.3 us Jsu]) and 



v' = 24 4- 30 [3l| have been reported, the spectroscopic constant u e is known [53|, |54| , etc. 

As for the 0+ state, a reference potential of three components (OM + RM + RM) has 
been constructed, but this time all parameters have been varied, except the equilibrium 



separation R e = 2.3893 A 47], which was kept fixed. The parameters of the PEC and the 
transition moment can be fitted independently, since the radiative lifetimes do not depend 
on the position of the levels. In its essence, the fitting procedure was the same as described 
in the previous subsection. As can be seen from Table 3 and Figs. 3 and 5, a good fit with 
the experimental data has been achieved for the l u excimer state as well. The transition 
moment curve shown in the inset of Fig. 2 is similar to that reported by Madej and Stoicheff 



3l|, but it falls to zero (l u —>■ 0;j~ transition is forbidden in the separated-atom limit) more 



slowly as r — > oo. 



III. FRANCK-CONDON FACTORS AND RADIATIVE LIFETIMES 

Calculation of Franck-Condon factors for the fixed pair of PECs and known transition mo- 
ment is a routine but rather time consuming task that has to be performed very accurately. 
In the present case, even more computational work is required because this demanding pro- 
cedure is used for fitting purposes. Fortunately, there are some possibilities to reduce the 
amount of computations. First, since the ground state reference potential is fixed, one has 
to calculate the related wave functions only once. Second, for any reference PEC under 
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examination, one has to solve the Schrodinger equation only at predefined points, e.g., at 
the abscissas of the relevant Gaussian quadrature formula. 

Thus, 640 wave functions for the ground state reference PEC have been calculated in 
the range E e [0,4.7 eV] with a variable energy step from 0.001 to 20 meV. Few examples 
of these eigenfunctions, all having the asymptotic form \l/(r) m 2cos((^3 — kr) according to 
Eq. (13), can be seen in Fig. 6. Of special interest might be the top graph in this figure, 
because it illustrates some fundamental findings of the quantum scattering theory at very low 



energies, i.e., where E = Ck 2 — > 0. Since tan ip 3 ~ — (fcr ) _1 [58J (r is the scattering length), 
and, according to Levinson theorem {6oj], the phase shift 5(k) ~ mr (n = 7 for the PEC 
in question) as k — > 0, one comes to the following expression for the wave function in this 
region: \?(r) ~ 2k(r — r ). This is explicitly demonstrated by a dotted line in the top graph 
of Fig. 6, where the vertical dotted line indicates the position of the characteristic parameter 
r . Naturally, such a simple linear coordinate dependence appears only at distances where 
the potential well becomes insignificant, i.e., V(r) ~ 0, but, on the other hand, r -C 

The scattering wave functions for the ground state reference PEC, along with the full set 
of bound state wave functions (only 7 in total), have been used as the basis to accurately 
ascertain all details of the Franck-Condon factors. Throughout the whole domain, the wave 
functions have been calculated at the abscissas of the 5-point Gaussian quadrature formula 
related to the intervals of 0.02 A width. Some results of these calculations are shown in 
Figs. 4, 5, and 7 to 9. In Fig. 7 one can see the Franck-Condon factors for bound-free 
transitions from the selected vibrational levels of both 0+ and l u reference PECs, while the 
details of the Franck-Condon spectrum from the highest level (v' = 29) of the 0^ state are 
demonstrated in Fig. 8. Note that the actual calculated probability density distributions 
are shown on both figures. From these data one can easily ascertain the probability of 



spontaneous emission (p v >) as well as the radiative lifetime 



t v i = 1/pv') of the levels. As 



explained in handbooks on quantum mechanics (see, e.g., 581 ]) 

oo 

4 



TV 



^{J2(E V , - E v „f[ J ^Ar)^ vl {r)dr} 2 (14) 



+ J {E v ,-Ef[ J V(E,r) f i(r)$ v ,(r)dr} 2 dE}, 
o o 

where E v > and $> v '(r) denote the initial (fixed) energy level and its wave function, while E v i. 
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and E are the discrete and continuous energy eigenvalues, respectively, with corresponding 
eigenfunctions ty v »(r) and ty(E,r). The calculations can be conveniently carried out using 
atomic units, i.e., taking H — 1, the velocity of light c = 137.03604, and measuring energy in 
Tartrees (1 Hartree = 27.2116 eV). Conversion to the SI frequency unit (Hz) is elementary 



59|: 1 Hz = 24.188843-10" 18 a.u. 



What are actually depicted in Figs. 7 and 8 are the energy integrands of the second 
term of Eq. (14) for the selected levels (including the factors and (E v > — E) 3 ). The 
low-energy part of the calculated spectra for the 0+ state (see the top graphs in Figs. 7 and 
8) nicely agrees with the experimental results by Moller et al. [24j], which confirms their 
assignment to the inner turning point region of the excimer's high-lying levels. On the other 
hand, this is an evidence of the validity of the ab initio PECs by Fernandez and Koch 4^]. 
As demonstrated in Fig. 7, the maximum of the second continuum is expected near 9.8 eV, 
which is also in full agreement with experimental observations. 

Another interesting result can be seen in the bottom graph of Fig. 8, which is related 
to the region where E — > 0, and, consequently, a very small energy step (0.001 meV) has 
been used. Namely, since ty(E,r) — > as E — > 0, the high-energy part of the probability 
density distribution looks like cut off. Indeed, the squared wave function of the v ' = 29 level, 
as needed, has 30 maxima, while only 28 maxima are seen for the related Franck-Condon 
spectrum in Fig. 8. Naturally, this simply means that one has to take the bound-bound 
transitions also into consideration. These contributions for the selected vibrational levels of 
both 0+ and l u excimer states are shown in Fig. 9. The discrete sets of 7 points (v" = 0-^6) 
there may seem to be located somewhat irregularly, but in fact, their positions are by no 
means accidental. On the contrary, as demonstrated in Figs. 4 and 5, the total integrated 
transition probability curves according to Eq. (14) are nice and smooth, although their 
bound-bound and bound-free constituents as if show some roughness. Such a behavior of the 
total probability is, of course, not only expected but even required. Indeed, this is nothing 
else but an indication of the completeness of the basis and the correctness of normalization 
of the wave functions. The overall transition probability for all levels should be exactly the 
same, if one puts fi(r) = 1 and ignores all factors in Eq. (14). Smoothly changing total 
probabilities simply reflect the smooth coordinate dependencies of the transition moment. 
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IV. CONCLUSION 



For any diatomic system, one can construct an exactly solvable multi-component reference 
potential based on the available experimental data. In this paper, we described a possible 
strategy to achieve this goal, which is analytically simple and computationally straightfor- 
ward. It is based on composing the reference PECs of several smoothly joined Morse-type 
potentials, and this was not an accidental choice. One might think that there are lots of 
alternatives, but this is not quite the case. Indeed, for any shape invariant potential one can 



easily find the exact solution of the related Schrodinger equation [44|. The point is, however, 
that they are already the physically correct linear combinations of special solutions, which 
can be easily ascertained only because the specific analytic form of the potential remains 
the same in the whole physical domain. The situation changes dramatically, if the potential 
consists of several analytically different pieces. Then it is often possible to easily ascertain 
one special solution but rarely both of them. A well-known exception is the piece-wise lin- 
ear potential possessing two linearly independent solutions in terms of Airy functions 61], 
nowadays available as standard functions in math-oriented programming languages. A use- 
ful method of solving the Schrodinger equation for a piece-wise linear reference PEC has 
been worked out long ago by Gordon [62 1. 

Another well-known example for which the two linearly independent solutions can be 
given in a simple analytic form is the Morse-type PEC analyzed in this paper. The classical 
Morse potential 43j itself is a reasonable approximation, but by introducing several com- 
ponents one can get a much more realistic description. Differently from a piece-wise linear 
potential, the components can be smoothly joined, preserving the continuity of the PEC and 
its derivative. The two linearly independent solutions of the related Schrodinger equations 
can be always given in terms of the well-studied confluent hypergeometric functions, and 
their correct linear combinations are determined by the boundary conditions and continuity 
requirements. 

The result of the fitting procedure depends not only on the theoretical methods and 
computational techniques applied, but also on the reliability of the experimental data and 
the constraints adopted. For example, following Fernandez and Koch [4^], we fixed the 
equilibrium nuclear separation R e = 2.3893 A for both 0+ and l u excimer states, but we 
cannot claim this to be the conclusive value for R e (slightly different results have been 
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reported by other authors). Since the transitions from the bottom of the excimer states fall 
into the repulsive wall of the ground state, the position of the second emission continuum 
strongly depends on R e . It means that if one assumes a different value for this parameter, 
the ground state PEC has to be changed as well. Another issue which probably needs further 
confirmation is the absolute numbering of the observed vibrational levels 30J, because the 
resulting PECs are rather sensitive to their changing. 

Nevertheless, as we demonstrated, a good fit with the experimental data has been 
achieved, and this is a strong argument for the reliability of the PECs obtained in this 
paper. The potentials involve a wide range of nuclear separations (r > 1.9 A) and they 
can be used to study the most important spectroscopic features of Ar^ excimers, including 
the details of their emission continua. In addition, they are expected to be useful for the 
analysis of the relaxation dynamics and the time-resolved emission spectra of the excimers. 
However, let us recall once again that these approximate PECs and the relevant transition 
moments have been deduced on the relatively modest basis of the available experimental 
data, following the relatively "soft" criteria stated in Section II. Hopefully, the results of 
this work can stimulate further experimental research to reveal much more details about the 
properties of the RG excimers. 
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Figure captions 



Fig. 1. Four-component reference potential (solid line) for the ground electronic state of Ar2 
in comparison with two ab initio potentials. The open circles in both graphs (and in 
;he inset) correspond to |47j, while the dotted curve in the upper graph is taken from 
481 ]. The same reference PEC (solid line) is depicted in both graphs, but essentially 
different energy scales are used for them. All components have the well-known analytic 
form of the Morse potential, but the ordinary Morse approximation is used only in 
the central range r e [X 2 , X 3 ] (see the lower graph). Two pseudo- Morse components 
for the regions r < X\ and r e [X!,X 2 ], respectively, have been introduced, while 
the long-distance range r > X 3 is approximated by a reversed Morse potential. The 
inset demonstrates that the RM approximation (dotted curve) does not contradict to 
the inverse power series expansion — Cg/r 6 — Cg/r 8 — Cio/r 10 (dashed curve). The 



parameters C6,C§,and C\q are taken from [171 ]. while the RM parameters for this 
specific fit are as follows: R max = 12.45 A, V = -0.1268 meV, D = -1.799E-4 meV, 
and a = 0.66 A^ 1 (V(R) = V + D * [exp(-a * (R - R rn )) - if). 

Fig. 2. Reference potentials for the three electronic states of Ar 2 examined in this paper. The 
parameters of the ground-state potential in the lower graph (the same as in Fig. 1) 
have been determined from the least squares fit to the ab initio potential by Fernandez 
and Koch [47], while the curves for the excimer states 0+ and l u have been constructed 
with the help of the fitting procedure described in Sections II B and C, respectively. 
The corresponding l u — > 0^" and 0^ — > 0+ transition moments are also shown in the 
upper graph (note that different scales are used for them), where the dotted line has 



been taken from 47J] 



Fig. 3. A more detailed depiction of the three-component reference potential for the 0^ state 
(upper graph), and the calculated vibrational energies for both excimer states in com- 
parison with the observed data (lower graph). As in the experiments, the position of 
the levels is measured relative to the zeroth level of the ground state. The v' — 31 
,ve. of the U state wa 8 taken f r0m 0, al. „t h e r expentnenta. data ate those front 



30|. 



Fig. 4. Upper graph: calculated probabilities of bound-free and bound-bound transitions 
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starting from the vibrational levels of the 0+ state. Lower graph: corresponding ra- 
diative lifetimes of all levels. The lifetime of the zeroth level (ro = 4.19 ns) practically 



coincides with the experimental result (r = 4.20± 0.13 ns) from 56|. 



Fig. 5. The same as in Fig. 4 but for the l u excimer state. Again, the lifetime of the zeroth 
level (tq = 3.15 fxs) is close to the experimental one (ro = 3.2± 0.3 fis) from j^ . 



Good agreement between the theoretical and experimental 3jj lifetimes for the levels 
v' = 24 30 is also demonstrated (see Table 3). 

Fig. 6. A selection of the scattering (E > 0) wave functions for the ground state (for correct 
energy normalization they have to be multiplied by F = (^ny/EC^j , where C = 
The upper graph corresponds exactly to the top of the tiny articial potential 
barrier (see the explanations in Section II B and Table 1), i.e., E = 1.029793-10~ 7 eV. 
Due to the extremely small E = Ck 2 , one can see nearly linear coordinate dependence, 
\I/(r) ~ 2k(r — ro) (ro is the scattering length), in the range where V(r) « but 
r ^ ^ ne three lower graphs correspond to E = 1 meV, E = 100 meV, and E = 4 
eV, respectively, when the wave function more or less rapidly achieves its asymptotic 
free- wave form. 

Fig. 7. Demonstration of the calculated Franck-Condon factors for the bound-free transitions 
from the selected vibrational levels of the excimer states. One can infer that the first 
and the second emission continua (with maximum near 9.8 eV) are formed just where 
expected, while the low-energy part of the spectrum nicely agrees with the experimen- 



tal results by Moller et al. 



241 ] under selective synchrotron radiation excitation. 



Fig. 8. Demonstration of the details of the Franck-Condon spectrum for the bound-free tran- 
sitions from the highest level of the 0+ state (shown also in Fig. 7). The lower graphs 
begin exactly at the energies where the upper graphs end. The bottom graph (where 
open circles mark the actually calculated intensities) is cut off at E = 0, and therefore 
the overall spectrum shows only 27 zeros (instead 29). The "missing" part corresponds 
to bound-bound transitions. 

Fig. 9. Depiction of the Franck-Condon factors for bound-bound transitions from the selected 
vibrational levels of the 0+ (upper graph) and l u (lower graph) excimer states of Ar2- 
The related atomic levels 3 Pi and 3 P2 are shown by the arrows. Location of the points 
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corresponding to the highest levels may seem irregular, but in fact, their positions are 
determined by rigorous quantum mechanical sum rules (see also Figs. 4 and 5). 
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TABLE I: Parameters of the reference potential for the ground electronic state of Ar 2 . The positions 
of the vibrational levels (in meV, relative to the bottom of the potential well) are also given. 
k Type V k (meV) D k (meV) a k (1/A) R k (A) Range (A) 

PM -67.83414 0.0763369 1.708238 5.224252 r < 2.9105 

1 PM -15.42556 0.137657 2.293931 4.443077 r G [2.9105, 3.35] 

2 OM -12.0866 11.85914 1.729656 3.7769 r £ [3.35, 4.4891] 

3 RM 1.029793E-4 -1.029793E-4 0.6877276 12.4891 r > 4.4891 

v" 1 2 3 4 5 6 

This work 1.8485 5.0741 7.6651 9.6403 10.9256 11.7307 12.0699 

Ref. 46 1.8127 4.9470 7.4416 9.3224 10.6342 11.4488 11.8691 

Ref. 30 1.8263 5.0102 7.5619 9.4923 10.8462 11.6918 
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TABLE II: Parameters of the reference potential for the 0+ state of Ar2. The vibrational energies 
given in meV are measured relative to the bottom of the potential well (this work), while the 



values in cm 1 are transition energies relative to the zeroth level of the ground state, as in 
The calculated spectroscopic constants are as follows: io e = 296.26 cm -1 , D e = 6128.3 cm -1 . 
k Type V k (meV) D k (meV) a k (1/A) R k (A) Range (A) 

1 OM 10864.76 1166.681 1.66221 2.3893 r < 3.0407 

2 RM 11624.65 -0.6649193 1.654645 4.862433 r e [3.0407, 4.16] 

3 RM 11623.59 -3.021853E-9 3.588833 7.00 r > 4.16 

v' 01234567 
E v ,(meV) 18.2936 54.4470 90.0222 125.0192 159.4380 193.2786 226.5409 259.2250 

v' 8 9 10 11 12 13 14 15 

E v ,(meV) 291.3308 322.8584 353.8078 384.1790 413.9720 443.1867 471.7515 499.7937 

v' 16 17 18 19 20 21 22 

E v ,(meV) 527.2506 553.8589 579.4339 603.7447 626.7231 648.2718 668.2571 

^(cm- 1 ) 92766.9 92940.7 93101.9 

Exp., Ref. 30 (cm- 1 ) 92769.3 92935.9 93093.5 

v 1 23 24 25 26 27 28 29 

E v r(meV) 686.6108 703.2166 717.9872 730.8195 741.5951 750.1805 756.3843 

E„-(cm^) 93249.9 93383.9 93503.0 93606.5 93693.4 93762.6 93812.7 

Exp., Ref. 30 (cm" 1 ) 93241.2 93377.6 93501.6 93610.8 93701.3 
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TABLE III: Parameters of the reference potential for the l u state of Ar2. The vibrational energies 
given in meV are measured relative to the bottom of the potential well, while the values given 
in cm -1 are the transition energies relative to the zeroth level of the ground state. The corre- 
sponding experimental values are from 



301 ] . except E31, which was taken from 57] . The following 
spectroscopic constants were obtained for the \ u state: uj e = 287.30 cm -1 and D e = 5929.6 cm" 1 . 
k Type V k (meV) D k (meV) a k (1/A) R k (A) Range (A) 

1 OM 10813.17 850.0 1.888481 2.3893 r < 3.39001 

2 RM 11548.17 -0.005519566 1.654645 6.390011 r £ [3.39001, 5.5] 

3 RM 11548.35 -1.515516E-10 0.730545 20.0 r > 5.5 



v' 





1 


2 


3 


4 


5 


6 


E v >(meV) 


17.7170 


52.5912 


86.7190 


120.1005 


152.7356 


184.6243 


215.7667 


v< 


7 


8 


9 


10 


11 


12 


13 


E v >(meV) 


246.1627 


275.8123 


304.7156 


332.8725 


360.2830 


386.9472 


412.8650 


v' 


14 


15 


16 


17 


18 


19 


20 


E v ,(meV) 


438.0364 


462.4614 


486.1401 


509.0725 


531.2584 


552.6980 


573.3912 


v' 


21 


22 


23 


24 


25 


26 


27 


E v/ {meV) 


593.3381 


612.4904 


630.8586 


648.2035 


664.3098 


679.0578 


692.3569 


E v / (cm.^ 1 ) 






92384.2 


92524.1 


92654.0 


92772.9 


92880.2 


Exp., Ref. 30 (cm- 1 ) 






92386.7 


92524.5 


92653.1 


92771.4 


92879.1 


Rad. lifetime (ns) 








182.1 


169.0 


159.3 


152.9 


Exp., Ref. 31 (ns) 








173 ± 17 162 ± 14 167.2 ± 9.0 161.7 ± 9.4 


v' 


28 


29 


30 


31 


32 


33 




E v ,(meV) 


704.0696 


714.0975 


722.3222 


728.6210 


732.8645 


734.8840 




E v / (cm" 1 ) 


92974.7 


93055.5 


93121.9 


93172.7 


93206.9 


93223.2 




Experiment (cm" 1 ) 


92974.8 


93056.9 


93123.7 


93171.0 








Rad. lifetime (ns) 


149.6 


149.4 


152.1 


158.0 


167.5 


185.6 





Exp., Ref. 31 (ns) 157.4 ± 9.6 155.7 ± 9.0 166 ± 11 
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